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We analyzed the Josephson bifurcation amplifier (JBA) readout process of a superconducting 
qubit quantum mechanically. This was achieved by employing numerical analyses of the dynamics 
of the density operator of a driven nonlinear oscillator and a qubit coupled system during the 
measurement process. In purely quantum cases, the wavefunction of the JBA is trapped in a 
quasienergy-state, and bifurcation is impossible. Introducing decoherence enables us to reproduce 
the bifurcation with a finite hysteresis. Moreover, we discuss in detail the dynamics involved when 
a qubit is initially in a superposition state. We have observed the qubit-probe (JBA) entangled 
state and it is divided into two separable states at the moment of the JBA transition begins. This 
corresponds to "projection". To readout the measurement result, however, we must wait until the 
two JBA states are macroscopically well separated. The waiting time is determined by the strength 
of the decoherence in the JBA. 

PACS numbers: 85.25.Cp, 05.45.-a, 85.25.Am, 03.65.Yz,42.50.Lc, 
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The readout of superconducting qubit states with the 
Josephson bifurcation amphfier (JBA) technique pro- 
vides non-destructive and high visibihty readout. There- 
fore, now it is widely and successfully used in actual ex- 
periments [l|. Mathematically, a JBA is described as a 
driven Duffing oscillator It enhances a small differ- 
ence in operation conditions by utilizing the bifurcation 
phenomenon. Under an appropriate driving force, a clas- 
sical nonlinear oscillator becomes bistable [2] . One stable 
state has a small amplitude (low-amplitude state), and 
the other has a larger amplitude (high-amplitude state). 
The critical driving force /c or the critical detuning 5c for 
the transition between these two states is very sensitive 
to small changes in the operational parameters of the os- 
cillator. For example, when we increase or decrease the 
driving force continuously, the amplitude of the oscilla- 
tion behaves hysteretically as shown in Fig. [T] When 
using a JBA as a qubit state readout probe, the JBA de- 
tects a small change depending on the qubit state. How- 
ever, the quantum-mechanical behavior of the JBA read- 
out process has not been established theoretically. This 
is because the bifurcation phenomenon can be discussed 
only for classical oscillators, and is impossible from the 
view point of pure quantum mechanics for an isolated 
system Q. A classical analysis gives no information on 
entanglement between the qubit and the probe (JBA) or 
the decoherence in the composite system, although all 
the quantum properties (projection, measurement back- 
action, etc.) in the readout are contained in such infor- 
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FIG. 1: Hysteretic behavior of the oscillation amplitude 
a{t) of a JBA as a function of the driving amplitude f{t) 
(schematic). 



mation. A quantum mechanical analysis is indispensable 
if we are to understand the readout process. 

In this letter, we analyze the quantum-mechanical time 
evolution of a JBA, and clarify how a bifurcation ap- 
pears in an actual situation. Moreover, we invesitgate 
what happens during the process of the qubit state read- 
out with a JBA by analyzing dynamics of the qubit-JBA 
composite system. 

In a highly quantum-mechanical JBA case, tunneling 
between classically stable states destroys the criticality in 
a classical oscillator. This type of phenomenon has been 
precisely discussed in ^ and in references therein, that 
the charging energy of the JBA ( ~ 2e^/C, where C is the 
effective capacitance in the JBA circuit ) is comparable to 
the energy barrier (^ the nonlinearity introduced below) 
between two stable states, and decoherecne is negligi- 
bly small. Howevr, actual JBA measurements are made 
with more classical conditions. Rigo et al. [3] investi- 
gated such an oscillator with a semi-classical trajectory 
analysis. In order to obtain quantum information more 
directly, here, we analyze the time evolution of a JBA 
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and a qubit during the readout process. 

A JBA can be modeled as an anhormonic oscillator in 
a rotating frame approximation with a Hamiltonian; 

Ht ^ {n^uj)na + ana'^ -^f{a^ +a) (1) 

where, at(a) is the creation (annihilation) operator of 
the Josephson plasma oscillation. Ua = a'' a, and fl is 
the linear resonant frequency of the JBA oscillator, lu 
is the driving frequency, which is slightly smaller than f2 
by the detuning S = — to. / is the driving strength, 
and a(> 0) is the nonlinearity. In a classical approxima- 
tion, this model shows the bifurcation in an appropriate 
parameter region. However, for a quantum-mechanical 
junction with [a,a^] = 1, the transition from |G)j (low- 
amplitude state) to (high-amplitude state) or, from 
\E)} to \G)j is impossible. 

The quasienergy-states (eigenstates of the Hamiltonian 
Eq. ([T]) in the rotating approximation) are easily calcu- 
lated and it is found that eigenstates never cross when the 
driving strength / is changed adiabatically. This means 
that if the JBA is initially in the ground state without 
driving, it never moves to the high-amplitude resonant 
state even if we increase the driving field because the 
JBA state only moves along the initial quasienergy-state 
and never jumps to the quasienergy-state which the high- 
amplitude resonant state belongs. 

Therefore, we expect that when a transition between 
quasienergy-states is caused by perturbation from out- 
side the system the bifurcation phenomenon is repro- 
duced. This is the case when decoherence is introduced 
into the present model. Here, we only take into account 
of the decoherence caused by a bath coupled to the JBA 
because decoherence that directly attacks the qubit is 
not limited to the readout process. Even for this model, 
indirect decoherence via the JBA occurs in the qubit. 

For example, we introduce linear loss in the oscillator 
(JBA) . The time evolution of the system (qubit- JBA) is 
governed by a Liouville equation: 

^^-[p^H] + ^{2apa^ -a^ap^pa^a), (2) 
at I 2 

where p is the density operator of the system, and T is 
the relaxation rate due to the linear loss in the JBA. The 
Q- value is given by fl/T. 

First, we show a numerical example of JBA dynamics 
without a qubit in Figs. [2]and[3l Here, the parameters 
used are S = 0.00717, a = 8 x lO'^n, Q = 2500, and / 
is operated as ^ 0.02512 — > 0. These parameters are 
similar to those used in actual experiments However, 
6 and F are a factor of 10~^ times smaller than real cases 
in order to emphasize quantumness and discuss the influ- 
ences of decoherence. Even if we use different parameters 
we obtain qualitatively same behaviors for a JBA with 
similar 6/T ratio value. 

Figure [31[a) approximately corresponds to the square 
of the JBA amplitude shown in Fig. [T] So, we can see 
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FIG. 2: (Color online) The trajectory of the quantum ex- 
pectation value ((^-j-^), (^"27^)) of the oscillator (complex) 
amplitude when the driving force / is operated as shown in 
Fig. [1] The starting point is \G}j without driving, and the 
right convergence point is \E)j. 
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FIG. 3: (Color online) Time variations of some quantum ex- 
pectation values when the driving force / is operated as shown 
in Fig. [1] (a) Number of bosons excited in JBA. (b) Purity 
of the JBA. (c) Fluctuation in JBA amplitude. 



that our calculation with decoherence reproduces the bi- 
furcation phenomenon well. We find that the critical 
driving /c is approximately O.OIIfi. Once the driving ex- 
ceeds this /c, the behaviors of the JBA is the same not 
depending on the maximum (/ = 0.025il for Fig. [3|) driv- 
ing strength. Moreover, our calculation provides a lot of 
quantum-mechanical information about the JBA transi- 
tion. Figurc[3jb) shows the time variation of the purity of 
the JBA state. Increasing the driving force /, we found 
that the purity decreases abruptly (t — t^). This cor- 
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responds to the beginning of the transition from \G)j 
to \E)j of the JBA. This is a manifestation of the fact 
that the transition needs an intense emission/absorption 
of energy to/from an external energy bath. This energy 
transfer is incoherent. After the rapid decrease, the pu- 
rity recovers to some extent and the JBA approaches to 
the classically stable state \E)j {t = t^). Since \E)j 
is a meta-stable state (stationary point of the classical 
Hamiltonian ), dragging JBA into the state by decoher- 
ence (linear loss) leads to the recovery of the purity. How- 
ever, the purity does not reach unity because it is not a 
true ground state. The fluctuation in the JBA amplitude 
is plotted in Fig. [3Kc). We can see a divergence of the 
fluctuation at the moment of the rapid decrease in the 
purity {t = Th). This suggests that this JBA transition 
between \G)j and \E)j is one of a phase transitions in 
bosonic systems with many degrees of freedom. 

Now we discuss the criterion of decoherence that de- 
termines whether a bifurcation is observed or not. From 
the above analyses we know that there is no critical value. 
When the decoherence is very small (F < S), the speed 
of the transfer from |G)j to |£')j becomes exponentially 
slower as (schematically) exp[— ?7(5/F], where 77 is a nu- 
merical factor of the order of unity. 

Information about the qubit state is transferred to the 
probe (JBA) through the formation of an entanglement 
between the qubit and the probe. What we actually ob- 
serve is the macroscopic state of the JBA, and merely 
postulate the qubit state. Therefore, the process by which 
the entanglement is formed and split into separable states 
due to decoherence ("projection") is very important for 
understanding the readout process Q . 

The qubit- JBA composite system is approximately ex- 
pressed by the Hamiltonian 



-(e(7z + ^cr^ 



(3) 



where Hq is the Pauli operator representation of the 
qubit. k is the interaction constant between the qubit 
and the JBA. The qubit state (a^) slightly changes the 
effective detuning 5 + ka^, resulting in a change in the 
critical value fc- By detecting the change in /c, we can 
distinguish the qubit state, i.e., whether is 1 or -1. For 
a flux qubit, the eigenstates of ctz are the two flux states, 
e is the bias provided by an external applied magnetic 
fleld, and A corresponds to the tunneling energy between 
two flux states. 

For the qubit-JBA coupled system, we carried out cal- 
culations similar to those without a qubit shown above. 
The qubit readout process is well understood by employ- 
ing knowledge of the quantum behavior in the time evo- 
lution of the JBA without a qubit that we have already 
discussed. 

We show a numerical example of the dynamics during 
the qubit readout process in Fig. 01 JBA parameters 
are the same as for the above example. The initial state 
is a separable state; {-^\g)q + -^\e)q) (g) \G)j, that is, 

the qubit is in a superposition. Here, \g)q and |e)q are 



the ground and excited states of the qubit, respectively. 
Qubit parameters are e = 0.2O, A/e — 1/2. The coupling 
between the qubit and the JBA is set at k — O.OOlil. The 
driving force / is increased from to 0.01 2il ( slightly 
larger than fc of the JBA) and maintained. This pa- 
rameter set gives a typical behavior of successful qubit 
readout. 

The Q-representations of the JBA state Trq[/5] are 
shown in Fig. |4l where p is the density operator of the 
qubit-JBA coupled system, and Trq[- ■ •] denotes taking 
partial trace about qubit degrees of freedom. In the read- 
out we can distinguish two peaks appearing in Fig. [IJd), 
which is the final stage of the readout. These peaks con- 
stitute an incoherent mixture, so they correspond to two 
possibilities in the measurement result. 

To discuss the entanglement between the JBA and the 



qubit, we adopt E = Tr _ Xr 



{TrqW 



as a mea- 



(Trq[p])^ 



sure of the entanglement. The reduction in Tr 

is the purity decrease in the reduced density operator of 
the JBA, that contains the decrease due to both decoher- 
ence and the entanglement formation. The reduction in 
Tr [p^] of the total system corresponds to the decrease 
due to decoherence. Therefore, E defined above shows 
the strength of entanglement. 

The time variation of the entanglement measure E is 
shown in Fig. [S] This process can be schematically ex- 
pressed as 



P(0) = \G)n{G\ ® + ^|e),) (^^^{g\ + ^^(. 

^ p(n) = i (|G),,|e)q + |G').,|5)q) (j(G|q(e| +,, (G'U(3|) 

P{r2) = i|G)j|e)qj(GU(e| + ^ |G') j |g)qj(G' 
^ p(r3) = i|G)j|e)q.j(G|q(e| + i|iJ),j|r;)q,,(iJU(gl. (4) 

Entanglement formation and "projection" correspond 
to the second {t = ti) and third {t = T2) fines of Eq. 
([4]), respectively. At this moment (T2), however, it is im- 
possible to obtain any information about the qubit from 
the observed probe (JBA) state because \G') closely re- 
sembles |G) in a classical mechanical sense (Fig. lU^b)) 
although quantum mechanically j(G'|G)j 0, namely, 
these two states are orthogonal. When we increase the 
driving force, one JBA state |G')j moves to \E)j. In con- 
trast, the other |G)j does not move significantly, (see. 
Figs. |ljc),(d)) Then (ra), we can easily distinguish \E)j 
or |G)j. This leads to a good postulation of the qubit 
state \g)q or |e)q, which brings us to the end of the read- 
out. 

The measure E is sufficiently quantitative for us to dis- 
cuss the time variation of the entanglement but it does 
not show the absolute strength of the entanglement. To 
estimate the absolute strength we can calculate the en- 
tanglement of formation for every eigenstate consisting 
the total system density operator p{t). For example, the 
time variation of the value of the most dominant eigen- 
state is quantitatively proportional to the behavior of 
E. However, it almost becomes unity when it reaches 
its maximum. The values for less dominant states also 
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FIG. 4: (Color online) Time-evolution of JBA during read- 
out. The figures show Q-representations of the JBA oscilla- 
tor states Trq[p] (in the rotating frame), (a) Beginning of 
the readout. The state of the total system is (schematically) 
p= |V'o)(i/'o|, with |i/>o) = l/V2(|ff),-H|e),)®|G)j. (b) Start- 
ing the transition. Entanglement formation and projection 
are carried out during this period, (c) During the transition. 
Entanglement has already been destroyed, (d) The entire sys- 
tem has become a mixture of classically correlated states. 



The backaction on the qubit caused by the measure- 
ment is induced as a result of the non-commutation 
relation between the qubit Hamiltonian and the inter- 
action Hamiltonian. When the qubit gap A is much 
smaller than other energies, the interaction commutes 
Hq. Therefore, the JBA readout causes only pure dephas- 
ing on the qubit. This does not pollute the measurement 
result because the measurement itself requires the pro- 
jection onto the basis. This is simply the condition of 
the "non-demolition measurement" . However, when A is 
not negligible compared with e, the measurement simul- 
taneously causes qubit relaxation. The non-commuting 
part induces coherent transition between \g)q and |e)q in 
the qubit. This coherent transition itself is not hurm- 
ful, but when such a transition is accompanied by deco- 
herence (linear loss), stochastic energy relaxation in the 
qubit accumulate and a finite error remains. In fact, in 
the numerical example shown above, the average (ct^) of 
the qubit deviates slightly (0.1%) from the initial value 
because of qubit relaxation. Stronger decoherence causes 
larger deformation in the readout result although it is of- 
ten much smaller than the deformation caused by other 
factors not discussed here, such as qubit relaxation as a 
reslut of decoherence directly attacking the qubit, even if 
we use 10^ times strong decoherence of JBA as in actual 
experiments. 



almost reaches unity. This means that the correlation 
between the qubit and JBA state becomes almost per- 
fect via the interaction between them. As a result an 
ideal JBA readout exhibits 100% visibility if the qubit 
relaxation discussed below is negligible. 
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FIG. 5: (Color on line) Time variation of the entanglement 
between the JBA and the qubit. Ti{i — 1,2,3) corresponds 
to those in Eq. |4l that is , ri: entanglement formation, T2: 
projection, T3: end of the readout. 



As described above, the state of the total system is al- 
ready divided into separable states |e)q|G')j and \g)q\G')j 
just after the transition |G")j — > \E)j starts. Therefore, 
"projection" itself is successful even if the transition takes 
much longer in the absence of sufficiently strong decoher- 
ence (here, linear loss F). However, we cannot distinguish 
|G)j and \E)j until the transition finishes. Then, the 
readout fails unless we can wait and maintain the JBA 
state until the transition is complete. 

In summary, we analyzed the quantum dynamics of 
the density operator of a system composed of a qubit 
and a JBA as the probe of the qubit state readout. From 
the analysis results, we have succeeded in extracting the 
essential feature of the JBA readout process of a super- 
conducting qubit. 
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